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I. 



INTRODUCTION 



Military forces have to be supplied with all kinds of 
items in order to complete their missions. Because many of 
the needed items cannot be purchased from outside suppliers 

at the time the forces need them, the items have to be 

/ 

purchased and stored earlier than actually needed. 

From the point of view of the ultimate user it would be 
best to store all such needed goods as near as possible to 
his own location, so that whenever an item is needed it can 
be delivered immediately. However, there are obviously 
constraints which make this course of action impossible. 

First, there is a budget constraint which allows only a 
certain amount of money to be spent on purchases of equipment, 
spare parts and other needed goods and, second, there is a 
space constraint because the supply support units are required 
to be mobile so that they can stay in close proximity to the 
combat units they support. In the case of direct support 
units for land forces the budget constraint is either non- 
existent or overwhelmed by the space constraint. 

In the Federal German Army these units are called BVP's 
and each consits of a fleet of trucks. They are equivalent 
to the Direct Support Units (DSU's) of the U. S. Army. The 
BVP's constitute the last echelon of the German Army's supply 
system and are not involved in buying goods from outside 
suppliers. Therefore the budget constraint on purchases is 
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nonexistent. However, the number of available trucks for 
a BVP is limited by its operational requirements and, to a 
certain extent, by the available budget for the armed 
services. Therefore any relevant inventory model has to take 
a space constraint into consideration when calculating optimal 
operating policies. 

Two models will be analyzed in this thesis. The first 
will be a multi-item deterministic lot size model subject to 
a space restriction which requires that the amount in 
inventory at any time cannot exceed the available storage 
space and a restriction that the time between orders for all 
items must be the same. The variables which are being 
optimized in this system are the time between orders for each 
item and the times between orders for different items. The 
notion of allowing orders for different items to arrive at 
different times will be referred to as time phasing of the 
orders for different items. The solution for this model is 

4 

obtained for the two-item case and the procedure for solving 
the n-item case is outlined. The second model is a multi-item 
deterministic order level system. It is subject to the same 
space and cycle time constraints as the lot size model. 
However, stockouts are allowed. The decision variables are 
the order levels of each of the items and the time phasing 
of the ordering of each item. This model is solved for the 
two-item case and discussed for the n-item case. 
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II. PREVIOUS WORK IN THE AREA 



A. CONSTRAINED LOT SIZE MODEL 

A multi-item lot size model with a single space constraint 
has been examined by Churchman [2] and Holt [6] . The model 
was defined in the usual form of deterministic inventory 
problems; that is, to minimize the total variable costs — 
reorder plus carrying costs — of the system. The space 
constraint required that the sum of the average space 
requirements for all items not exceed the total space available 
for these items. The assumptions which were made in the 
formulation of the problem were: 

1. a known constant demand rate; 

2. infinite replenishing rate; 

3. constant lead time; 

4. shortages are not allowed. 

The cycle lengths of all items were not required to be 
identical. In addition, the time phasing of each item was 
not addressed. 

A disadvantage of this model is that only the sum of the 
average inventories must be no greater than the space available. 
This requirement can be satisfied mathematically and yet not 
be met in reality because the possible different cycle lengths 
may cause the order for more than one item to arrive at the 
same time. 
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B. A CONSTRAINED ORDER LEVEL SYSTEM 



A multi-item deterministic order level model has been 
examined by Naddor [7] . This model used the first three 
assumptions of the above model by Holt and, in addition, 
required a fixed review interval for all items. Backorders 
were allowed however. The objective function for minimization 
was the total variable cost — carrying plus backorder costs — 
and the constraint required that the summation of the space 
requirements of the optimum order levels not exceed the total 
space available. 

The major difference between Naddor' s model and the one 
to be presented in this paper is that in his model all orders 
arrive at the same time, while the model being developed here 
jf^ermits time phasing of arrivals of orders of different items. 
His solution for the optimal order level for the n-item case 
is : 



c 0 . - g b . 
c - r, 2i ^o 1 

i *^i c. . + c- . 

li 2i 



where 

= order level of item i when lead time is zero; 
q^ = replenishment size of item i; 
c^ = the carrying cost of item i per unit of time; 
c 2i = t ^ le s h° rta 9 e cost of item i per unit of time; 
b^ = unit space required for item i; 

B = total warehouse space available; 
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and 



n c- . b . q . 

- o. - W t - B) 



K 2 rr 

d . q . 

l ^1 



n 
2 

. i c, . + c 0 . 
i=l li 2i 
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III. A CONSTRAINED MULTI-ITEM LOT SIZE SYSTEM 



The model which is developed in this chapter is an n item 
lot size model constrained by the amount of space which is 
available for storage of inventory. The following assumptions 
are made: 

1. Demand for item i is deterministic at a constant 
positive rate of r^ units per unit time. 

2. No shortages are allowed and the stock of an item will 
be immediately replenished whenever its inventory level 
reaches zero. 

3. The lead time is constant. 

4. The replenishment rate is infinite. 

5. One type of storage is usable by all items. 

6. The inventory carrying cost is proportional to the 
average amount in inventory at any time; its dimension 
is $ per unit-time. 

7. The order cost is a constant which is incurred each 
time an order is placed; its dimension is $. 

The variables which will be optimized to produce minimum 
cost are the length T of the reorder period for every item 
and the time x^ at which an order for the ith item arrives 
after arrival of an order of item number 1 (x^ = 0). The cost 
equation, as a function of the scheduling period for any 
single item under these assumptions, is given by 



C i (T) = c li r i T/2 + C 3i /T ; 



( 1 ) 



where c. . = the inventory carrying cost, and 



c 0 . = the unit reorder cost. 
3i 
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This notation is that of Ref. 7, page 58. For n items, the 
total cost can be expressed as 



n 



TC 





given that all items are using the same reorder period length. 
The space constraint can be formulated as 



b^ = the space occupied by a unit of the ith item; 
I^(t)= the inventory level of the ith item at time t. 

Since maximum inventory levels occur at the times at which 
orders are received, these times correspond to critical points 
relative to the constraint. Therefore the constraint needs 
only to be examined at these arrival times. This will result 
in n constraints for an n-item case because of the possibility 
that each item may arrive at a different time. 

To solve the problem let represent the ith Lagrange 
multiplier. Then the Lagrangian becomes 



n 

Z 

i=l 



b.I . (t) < B 

11 — 



for all t > 0 



(2) 



where 



B 



total space available 



L 



n 

Z 

i=l 



C ± (T) 



n 

+ Z 
i=l 




(3) 
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where C^(T) is given by equation (1) and represents the 
ith constraint function (describing (2) at each x^, 
i=l, 2 , . . . ,n) . The solution to equation (3) will be presented 
for cases which have two and n items. The two-item case is 
presented to illustrate the general solution procedure. 

A. SOLUTION OF THE TOO- ITEM CASE 

Equation (3) for the two-item case reduces to: 

L = c-^r^T/2 + c-^/T + c ;L2 r 2 T / 2 + c 22 /T + ^1 F 1 + ^2 F 2 * 

Typical fluctuations of the inventory levels of the items can 
be represented as shown in Figure 1. 

From Figure 2 it can be seen that the constraint function 
at time zero when item number 1 arrives will be 



and that the constraint function at time x 2 when item number 
2 arrives will be 



Conversion of inequalities into equalities for use in the 
Lagrangian yields: 



b l r l T + b 2 r 2 X 2 — B 



b 1 r 1 (T x 2 ) + t> 2 r 2 T £ B 




2 



B 



(4) 



F 2 = b l r l (T - x 2 > + b 2 r 2 T + V 2 



2 



B 



(5) 



where and V 2 are slack variables 



I(t) 




FIGURE 1: Fluctuation of inventory levels for two items for 

the deterministic lot size model. 
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FIGURE 2: Space occupied by two items when = b 2 r 2 
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Combining (3), (4) and (5) gives: 



L = c^r-jT/2 + c 31 /T + c 12 r 2 T/2 + c 32 /T 

+ X^(b 1 r^T + b 2 r 2 x 2 + V^ 2 - B) 

+ X 2 (b 1 r 1 (T - x 2 ) + b 2 r 2 T + V 2 2 - B) . (6) 



Taking partial derivatives of the Lagrangian with respect to 
V^, V 2 , x 2 , T, X^ and X 2 , setting them equal to zero, and 
simplifying, yields: 



Xi v i 



= 0 



X 2 V 2 



= 0 



^l b 2 r 2 *2 b l r l 0 



C ll r l //2 C 31 //T + C 12 r 2' /2 C 32 //T + X l b l r l 



+ ^ 2 b l r l + ^2 b 2 r 2 0 



(7) 

( 8 ) 

(9) 

(10) 



b l r l T + b 2 r 2 X 2 + 



V, 



- B = 0 



(ID 



b l r l(T-x 2 ) + b 2 r 2 T + 



V, 



- B = 0 



(12) 



Equations (7), (8) and (9) indicate that there are only 

two possible cases: 



17 



1 . X^ 7 ^ 0 and ^ 7 * 0 which implies that = V 3 = 0 ; 

2 . = X 2 = 0 , and the values can be arbitrary. 

When X^ j* 0 and X 2 ¥■ 0 , both constraints are binding. Solving 
equations (11) and (12) with = V 2 = 0 for T and x 2 
respectively, yields: 

T* = Btb^+b^J/tb^r^+b^r, + b 2 2 r 2 2 ) . (13) 

x 2 = Bb 2 r 2 /(b 1 2 r 1 2 + b 1 r 1 b 2 r 2 + • (14) 

Substituting T and x 2 into the total cost equation yields 
the minimum cost C where 



C o = B(c 11 r 1 + c 12 r 2 ) (bj^r^^ + b 2 r 2 ) /2D + D(c 3 l +c 32 )/A , 

and 

D = b l 2 r l 2 + b l r l b 2 r 2 + b 2 2 r 2 2 ; 

A = B(b^r^ + b 2 r 2 ^ • 



The ratio 



provides a convenient means of comparing ^ to T for this 

completely constrained case. If b^r^ = ' D 2 T 2 then 

* * 
x yi 

-| = J • If Vn < b 2 r 2 then 4 > j ■ If > b.,r 2 
* 

X 2 1 

then — * < y . 

T 

For the unconstrained case, equation (9) requires that 
if one of the X values equals zero, the other has to be zero 
too. The optimal unconstrained value of T, denoted as T , 
can be determined from equation (10) . It is 



T 



u 




C 32 ) 



C 12 r 2 



(15) 



Equation (15) is the same expression as derived by Naddor [7]. 
The value for is dependent, however, on the value of B. 

The value of B which corresponds to the dividing point between 
cases 1 and 2 is obtained by equating equations (13) and (15) 
and solving for B. This value is denoted as B^, and is given 
by equation (16) 



B 1 = 



D 

b l r l + b 2 r 2 



u 



(16) 



When B <_ B^ , equations (13) and (14) give the values of T* 

and x- . When B > B. , T = T and the value of x- is no 
2 1 u 2 

longer restricted to that given by (14) . 
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If B &2 where = ^^^4^2^) 



* 

then X 2 can be 



any value in the interval 0 £ X 2 _< T . This corresponds to 
B being so large that even if orders for both items arrive 
simultaneously they will not be faced with a storage 
restriction. When < B < B 2 , the storage constraint will 
not be restrictive provided that X 2 is so selected that the 
constraint is not exceeded at the time either item arrives. 

When the first item arrives (at t = 0) the total combined 
on-hand inventory consumes b^r^T^ + b2 r 2 x 2 amount of 
storage space. When the second item arrives (at t = X 2 ) the 
total combined on-hand inventory consumes space in the amount 
of b^r^(T u ~X 2 ) + ^ > 2 r 2^'u * T’he va ^ ue x 2 mus ^ selected 
so that the total combined on-hand inventory at either time 
does not exceed B. Thus if X 2 lies in the interval 



then both conditions will be satisfied when B^ < B <22* 

When B = B 2 the interval given by (17) simplifies to 
® 1 x 2 1 T u • As B decreases below B^r the upper bound 
decreases and the lower bound increases at the same rate. 

The interval reduces to a single value for x 2 at B = B^; 
the corresponding X 2 value is given by equation (14) . As a 
consequence of this behavior of the interval for B between 
B^ and B 2 , it follows that equation (14) gives an optimal 
value of X 2 for all B > B^ , and the only optimal value for 
B £ B^ 




(17) 
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B. SOLUTION OF THE n-ITEM CASE 



The Lagrangian for n items is: 
n 



n 



L = Z (c,.r.T/2 + c,./T) + Z A F. , (18) 

i=l 11 1 - 31 i=l 1 1 

with the constraint functions given below. These functions 

assume 0 < x~ < x 0 < • • • < x 

— z, z n 

b l r l T + b 2 r 2*2 + b 3 r 3 x 3 + ’ ' ' + ViA + v l = B 

+ b 2 r 2 T + b 3 r 3 (x 3‘ x 2 ) + ••• + VnW + V 2 = B 

b l r l^ T-x 3 ^ + b 2 r 2 (T - (X3-X2)) + b 3 r 3 T + ... + b n r n ( x n “ x 3 ) + v g 2 = B 



b 1 r 1 (T-x n _ 1 ) + b 2 r 2 (T (x^j-x^ ) + ... + b n r n (> Yi _X n-l ) + V n-1 ” B . 

b l r l'^V +b 2 r 2 (T -<V x 2> + ••• +b n r n T + V n 2 = B 



After forming the Lagrangian, taking partial derivatives with 
respect to x 2 , x 3 , . .., x n and V^, V 2 , ..., V n , and setting 
them equal to zero, it can be shown that either all of the 
Lagrange multipliers are zero or all slack variables are zero. 
The partial derivatives with respect to A^, X 2 , ...» A , 



21 



when the V^'s are all zero give back equality constraints at 



t = 0/X2».../X n . Solving these equations for T, X 2 , 
yields : 



* 

T 



B 



n 



E 

i=l 



b . r . 



1 1 



D 



• • • / 



X 



n 



/ 



( 19 ) 



* 



x. 

1 



B 



l 

E 

k=2 



Vk 



1 , 2 , 



» n 



( 20 ) 



where 



D 



n 

E 

i=l 



, 2 ; 

b. r . 

1 1 



+ 



n-1 

E 

i=l 



n 

b . r . z b .r . 
1 1 j=i+l ^ 3 



( 21 ) 



The resulting value of C Q , the minimum total variable cost, 
is : 



n 



n 



C Q = B ( Z c ii r (t>i r i+. ..+b n r n )/ D + D ( E c 31 /B (b 1 r 1 +...4b n r n ) . 



In the unconstrained case, optimal T is obtained from 




( 22 ) 



The n-item equivalents of B^ and B 2 are: 

D T 



B- 



u 



n 

E b . r . 
i=l 1 L 



( 23 ) 
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( 24 ) 



n 




And, if E >_ B 2 1 then T = and each x^ can take on any 

value in the interval 0 £ x i £ T u * When < B < then 
* 

the x^ ' s must be chosen in a way that the constraints are 
not violated. As B decreases below B 2 , the upper bound 
decreases and the lower bound increases. However, equation 
(20) has been obtained for the case where the constraints 
have to be met, i.e. where we had equalities in all of the 
constraint functions. Changing the equalities to inequalities 
will result in an infinite set of optimal feasible solutions. 
But x^ as given by equation (20) will always be in this set 
of optimal feasible solutions because the solution to the 
more restricted problem will always be an answer to the less 
restricted problem. 

The set of optimal solutions creates the ranges for the 
x^ values. They could be obtained using linear programming 
in the following way : 

1. Possible objective functions are 

a) min x. , i = 2,3,...,n, to get the smallest lower 
bound for x^ ; 

b) max x^, i = 2,3,...,n, to get the greatest upper 
bound for x^ . 

2. The constraints are all of those shown below 
equation (18) . 
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For solving problems which have n > 2 items the flow chart 
of Figure 3 is helpful and summarizes the preceding analysis. 

C. EXAMPLES 

Two examples are presented below to illustrate the model. 
Example 1 

Consider a warehouse which stores two different items 
whose parameters are: 

Item 1: Item 2: 

r^ = 200 units per year r 2 - 250 units per year 



If total available warehouse space is 400 cubic feet, what 

* 

are the optimal order period T for each item and the optimal 
* 

arrival time X 2 for item number 2 after an order of item 
number 1 arrives? 

Calculating and yields: 



b 




c 



^ = 5 ft^/unit 

= $ 2 per unit-year 



b 2 = 3 ft^/unit 

c 12 = $ 1 per unit-year 
c 3 2= $ 25 



= 0.41 years ; = 543 ft 



3 



Because B < B^ , both storage constraints are binding. 

Vt it 

Therefore we get T =0.30 years and x 2 = 0.13 years by 
using equations (13) and (14). 
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FIGURE 3 




Flowchart for solving lot size model problems 
having n > 2 items. 
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Example 2 



Suppose the data of example number 1 remains unchanged 
except for B. Assume B is now 600 ft . 

The values for and B^ are unchanged. But since B is now 
greater than B^ evaluation of B 2 is required. It is 
calculated to be 719 ft 3 . 

★ 

Because B-^ < B < B 2 r T = T^ = 0.41 years and x 2 lies 
in the interval given by (17); that is, 

0.12 years <_ x 2 <_ 0.25 years . 

Checking the feasibility at time t = 0 shows that, using the 

* 

lower bound as the x 2 value, the occupied warehouse space 

3 * 

would be 500 ft ; using the upper bound as the x 2 value, 

3 

600 ft would be needed. The respective space needs for 
time t = x 2 would be 600 ft 3 and 467 ft 3 . 

D. COMMENTS 

The constrained multi-item lot size system developed in 
this chapter differs from that of Holt [6] in the following 
ways : 

1. Holt's model allows a different cycle length for each 
item, while a common cycle length for all products is 
used in the above model. 

2. Holt's constraint required that the sum of average 
inventories be no greater than the space available, while 
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the constraint function used in the above model requires 
that the space occupied by the inventory at any time be 
no greater than the total space available. 

3. The time phasing between orders for different items is 
calculated in the above model, while the time at which 
an order arrives is not regulated by Holt. 

The second aspect is perhaps the most important difference 
because the actual space restriction in Holt's model is a 
nebulous notion. Realistically, the space constraint restricts 
the maximum on-hand inventory rather than the average level. 

It is interesting to note that the solution for the 
scheduling period can also be used to calculate a minimum 
amount of storage space needed given a specific scheduling 
period. The equation for the warehouse space in terms of the 
scheduling period is 

I 

n 

B = T D / l b . r • 
i=l 1 1 

where D is given by equation (21) . 
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IV. A CONSTRAINED MULTI-ITEM ORDER LEVEL SYSTEM 



A. THE MODEL 

The system which is analyzed in this chapter is an n-item 
order level or periodic review system which is constrained by 
the total amount of space which is available for storage of 
the items. The analysis of this system will be conducted 
under the following assumptions: 

1. Demand for item i is deterministic at a constant rate 
of r^ units per time. 

2. The period between replenishments or scheduling period 
is a fixed constant T. 

3. The replenishment quantity of the ith item raises the 
inventory at the time the quantity is received to the 
order level after backorders are made up. 

4. The replenishment rate is infinite. 

5. The lead time is constant. 

6. The same type of storage is used by all items. 

7. The inventory carrying cost is proportional to the 
average amount in inventory. Its dimension is $ per 
unit time. 

8. The backorder cost is proportional to the average shortage 
level. Its dimension is $ per unit time. 

9. The order costs are ignored because they are a fixed cost 
if the scheduling period is a prescribed constant. 

From these assumptions it can be seen that the lot sizes which 

are being considered are constant since 



= r . 



T 



/ 



where q^ represents the ith item lot size. 
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The variables which are to be optimized in this system 
are the order levels (S^ for the ith item) and the points in 
time at which each item arrives (x^ for the ith item, where 
x 1 = 0) . 

The cost for one scheduling period for the ith item as a 
function of the order level can be shown to be [7] : 



C(S i ) = 


c li S i 2 „ c 2i (r i' S i )2 
2r i T 2r i T ' 


where c, . = 

li 


inventory carrying cost. 


and c~. = 

2i 


inventory shortage cost (using Naddor’s 
notation) . 



Then the total cost equation for an n-item model can be 



expressed as: 




C (S^ , S 2 ! • 


2 2 
n c., . S . c 0 . (r . -S . ) 

O \ — V / 1 -f 2. 1 1 1 V 

• * ' n' . . 1 2r.T 2r . T 1 

i=l 1 1 



The storage constraint can be formulated as: 
n 

E b.I. (t) < B for all t = 0 , 
i=l 1 1 



where B 


= total space available, 


b. 

1 


= amount of space which is occupied by one unit 
of item i, 


I^t) 


= inventory level of the ith item at any point 
in time t. 
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As in the preceding model, the maximum amount in inventory 
of any item occurs immediately after the ordered item arrives. 
Therefore the arrival time is a critical time with regard to 
the space constraint. Because a constraint will be required 
for the time of arrival of each item, there will be n 
constraints for n items. 

Figure 4 illustrates the two possibilities in fluctuations 
of inventory for two items. Two cases are possible because, 
at the end of each scheduling period, an order for one item 
arrives and the other item is either on backorder or in stock. 
Care must be taken in formulating the constraints because when 
an item is on backorder its inventory level is negative but 
the amount of its occupied space is zero. 

B. SOLUTIONS FOR THE TWO- ITEM CASE 

The problem for the two- item case can be stated as 
follows : 



Minimize 



TC _ C 11 S 1 2 j C 21 (r l T ~ S l ) t C 12 S 2 + C 22 (r 2 T ~ S 2 ) 



2r x T 



2rjT 



2r 2 T 



2r 2 T 



(25) 



subject to the constraints 



b-^S^ + max[0, b 2 (S 2 - r 2 (T-x 2 ) ) ] <_ B ; (26) 

b 2 S 2 + max[0, b^(S^ “ r^x 2 ) J — B ( 27 ) 
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Case 1 



FIGURE 




Case 2 




: Fluctuation of inventory levels for two items for 

the deterministic order level model. 
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Constraint (26) reflects the combination of the constraint 
cases at t = 0. Whenever S 2 -r 2 (T-x 2 ) < 0 , which means that 

item number 2 is on backorder, this constraint reduces to the 
restriction that b-^S^ be less than or equal to B. Similarly, 
constraint (27) reflects the combination of cases at t = x 2 . 

1. Totally Unconstrained Case 

If the problem is unconstrained, then the optimal value 

for S-^ and S 2 can be obtained by simply taking the partial 

derivatives of equation (25) with respect to and S 2 , setting 

them equal to zero and solving for and S 2 . The solutions, 

denoted as S , and S are: 
ul u2 



S 



ul 



C 21 r l T 
C 11 + C 21 



(28) 



5 u2 



C 22 r 2 T 
C 12 + C 22 



(29) 



which are the same expressions as in Naddor's unconstrained * 
case [7]. Because B _> b i s u l + b 2 S u2 ~ B 3 ' even simultaneous 
arrivals are not space constrained and x 2 can assume any 
value in the interval 0 £ x 2 £ T . 

2 . Partially Unconstrained Case 

When B < B^ , then simultaneous arrivals would 
violate the space constraint. If B is not too restrictive, 
the problem could still be unrestricted relative to and/or 

S u2 if careful time phasing of the orders is done. 
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The individual constraints were: 



+ max [0 , ^2 < S 2 ~ r 2 ^ T-x 2 ^ 1 £ B 
b2S2 + max[0, b^fS^-r^^)] £ B 



Notice that S ^ and S 
space constraint if B 



u 2 will not necessarily violate the 

B 2 and X 2 is chosen carefully, where 



B 2 = max [b 1 S- ul , b 2 S u2 l 



(30) 



They will therefore remain optimal. 

If we assume that both b-^S^ < B and ^2 S u2 < B 
then the functions max[0, b 2 ( s u 2 “ r 2 ( T-x 2 ^ 1 and 
max[0, ( s u i ” r x x 2^1 at t = 0 and t = X 2 , respectively, 
could be positive. At t = 0 the value of the max function 
is not to exceed the total space available minus the space 
occupied by item number 1. A feasible positive max function 
at t = 0 satisfies: 



0 < b 2 (S u2 r 2^ T X 2 J) - B b l S ul ' 



or, in terms of X 2 , 



T — < x 0 < T 

r 2 2 " 



fu2 + B b l S ul 
r 2 b 2 r 2 



(31) 
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If x 2 < T - then b 2 (S ^2 ~ r (T-x 2 ) ) <_ 0 . If x 2 is 

greater than the right-hand bound then the total inventory 
at t = 0 would exceed B, the storage restriction value. 

When t = x 2 the value of the positive max function 
is not to exceed the total space available minus the space 
occupied by the optimum order level of item number 2. This 
means that 

0 •= b l< s ul' r l x 2> i B ' b 2 S u2 ' 



or, in terms of x 2 , 




B - b 2 S u2 
r l b l 



< x 2 




(32) 



In this case if 




then b i ( s u i " r i x 2 ^ 1 0 



and if x 2 < left-hand bound, the total on-hand inventory 
at t = x 2 would exceed the storage restriction. 

In comparing inequalities (31) and (32) it is clear 
that the space constraint will be only met if: 




B - b 0 S 0 
2 u2 



< x~ < T 



r l b l “ 2 - 





b.S . 
1 ul 



2 r 2 



(33) 



However, feasibility requires that the left-hand side of (33) 
be less than or equal to the right-hand side. 
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The two bounds of (33) are equal when B = B where 

X 



B x B 3 " B 0 ' 



( 34 ) 



and 



B_ = b,S , + b 0 S 0 ; 
3 1 ul 2 u2 



(35) 



b l r l b 2 r 2 

B. = [ r - — - — ~r~r — ~ — ] T . 

0 b l r l +b 2 r 2 



(36) 



Therefore, if B 2 £ B x £ B < B^ then the problem is 

* 

unconstrained with respect to and S 2 provided x 2 is chosen 
within the boundaries of inequality (33) . 

Simplifying (33) results in: 



Bo “ B Bo ~ B 

-4 < x 0 < T ^ 

b l r l - 2 - b 2 r 2 



(37) 



It follows that when B = B^, the range given by (37) is 

< 

0 £ x 2 £ T ; when B < B^ the range becomes narrower with the 
upper bound decreasing at a rate l/b 2 r 2 and the lower bound 
increasing at a rate -1/b^r^ per one unit decrease in the 
value of B. We note in passing that if B = b^S^ = £>2 S u2 ' 
the range (33) reduces to 



S , S 0 

-Hi < x, < T - -H 2 
r l - 2 - r 2 
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The condition resulting in B > max[b.S , ,b~S „1 is a case 

x 1 ul 2 u2 

where both constraints are binding for any B < B while the 

X 

max functions are both still positive for B near B^. 

3 . Constrained Cases 

a. B > B 0 and B A < B < B 
x 2 0 x 

If B^ > max[b 1 S ul ,b 2 S u2 ] and B < the 
constraints when the max functions are both positive are: 



b 1 S 1 + b 2 S 2 - b 2 r 2 T + b 2 r 2 x 2 < B ; 



b 2 S 2 + b 1 S 1 - b iri x 2 



< B . 



Subtracting the second inequality from the first gives: 



= [ 



b 2 r 2 



B, 



b l r l +b 2 r 2 



) T = 



b l r l 



( 38 ) 



Thus, x 2 is independent of B provided that both max functions 

* 

remain positive. Substitution of x 2 into either constraint 
gives : 



bfSi + b 2 S 2 B + t>^r^x 2 



< B + B 



0 



Therefore we have the very simple relationship that, whatever 
values of S-^ or S 2 are chosen, we must have 

( 39 ) 



h l S 1 + b 2 S 2 1 B + B 0 



36 



Now inequality (39) does not reflect the fact that the space 

constraints at t = 0 and t = x 2 would be violated if either 
* * 

bfSf > B or t>2 S 2 > B * Therefore, the following conditions 
must also be imposed: 

b l s l 1 B > (40) 

b 2 S 2 - B * (41) 

The form of the resulting constraint region is illustrated 
in Figure 5(a). The unconstrained optimal solution (S u ^,S u2 ) 
lies outside this region. Because the objective function is 
convex in and S 2 and the constraint region is convex, it 
follows that the constrained optimum lies on the boundary of 
the feasible region. 

If B is reduced to Bq we get the condition illustrated 

in Figure 5(b) and when B < Bq we get the result of Figure 

5(c). These cases will be analyzed in the next section. 

To search for the optimal solution when Bq < B < B^ 

we evaluate and S 2 using only one segment of the boundary 

* 

at a time. We begin by assuming that the constrained and 
★ 

S 2 lie on the line given by (39) with equality. We can then 
use (39) to solve for S 2 as a function of S^, substitute that 
function into the objective function, differentiate with 
respect to S^, set the derivative to zero and solve for S^. 

We get 
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FIGURE 5 : Feasible regions for the two-item deterministic 

★ ★ 

order level model where b-jS-j + b2$2 <_ B + Bg . 
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* ~ 

s i = s : 



C 21 “ C 22 [ b 2 ] + 



® + ^1 c 22 
°] t^] Ic r -] 
°2 b u2 



C 21 ^1 ^ c 22 

tc^] + [r^] [^] 

ul D 2 u2 



( 42 ) 



* ~ 

S 2 = s. 



B + B o b i - 

5 ‘V 1 • 21 



(43) 



If either or S 2 violates the b^S^ _< B or b 2 S 2 _< B 
constraints then the optimal constrained values lie on eithe 
of these latter constraint boundaries. 

A 

If is larger than B/b^ then the largest feasible 

value is B/b^. This then will be the optimal value of S^ 

* 

Substitution of = B/b^ into objective function and 

differentiating with respect to S 2 will give S^ 2 when the 
derivative is set to zero. Therefore, optimal constrained 
S 2 is the largest value the constraints will allow without 

ie 

reducing . This corresponds to the S 2 value at the inter 

* 

section of (39) and (40); that is S 2 = Bg/b 2 . 

★ 

Similar arguments lead to S.^ = Bg/b^ and 
* A 

S 2 = B/b 2 when S 2 is greater than B/b 2 . 

★ 

In summary, if B Q < B < B^ then x 2 is given by 

★ ★ 

equation (38). To determine S^ and S 2 : 

A A 

1. Calculate S^, S 2 . 

^ ^ ~ ~ 

2. If S 1 < B/b 1 , S 2 <_ B/b 2 then s i = s i ' S 2 = S 2 * 
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3. 


If 


A 

S 1 > 


B/b^ , then 


* 

S 1 " 


B/b 1 , 


S 2 B 0 //b l * 




4. 


If 


A 

S 2 > 


B/b 2 , then 


* 

S 1 = 


B 0 /bi , 


S* = B/b 2 . 






b. 


B 

X 


> B 2 and B 


1 B o 


• 










Figure 5(b) and 


5(c) 


emphasize 


that when B < 


B 0 


the 


constraint 


given by (35) 


lies 


outside 


the boundaries 


of 



the feasible region and is trivially satisfied by any solution 

in the feasible region. It follows from the arguments in the 

* * 

preceding section (a) that = B/b^ and S 2 = B/b 2 and 

we realize immediately that both max functions must be zero. 

The analysis of the preceding section showed that 



if 




/ 



then both max functions will remain zero. Substitution of 
* * 

and S 2 into this range gives 



B 



b l r l 



1 x 2 1 T 



B 

b 2 r 2 



(44) 



This range reduces to x 2 = B o //b l r l ( e< 5 uat i° n (38)) when 
B = Bq ; the range increases to 0 < x 2 < T as B is reduced 
to zero. 

Figure 6 summarizes the ranges of optimal x 2 for 
the case where B x > ®2 * Brom this figure it becomes obvious 
that x 2 = Eg/b^r^ °Ptimal regardless of the value of B. 
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FIGURE 6: The range of optimal as a function of B 

for the case where B x > B 2 . 



41 



c. 



B 2 > B x 



and B < B, 



In the constrained case where B 9 > B , further 

Z X 

reduction into cases where b,S , > b 0 S „ , b 0 S „ > b n S , 

1 ul 2 u2 2 u2 1 ul 

and b^S ul = b 2 S u2 are necessary. We will also make use of 



B]_ E min[b l S ul' b 2 S u2 ] 



(45) 



(1) B 2 bf^ui / b 2 B u2 ' B ^ B 2 • 

When B^ < B < B 2 then the first constraint 
to become binding is that at t = 0 . Its max function will be 
zero. An appreciation of what is taking place can be obtained 
by looking at the x 2 range given by (33) . When B 2 > B^ then 
the upper and lower bounds do not reduce to the same value 
(at B = B^) before the B - b l S ul term in the upper bound 
goes to zero. The resulting range for x 2 from (33) is now: 



b 2 S u2 

b l r l 



< x 2 < T 




(46) - 



* 

when = B/b^ is substituted for S^. The lower bound is 

only valid however as long as B B^ = b 2 S u2 * Tbe °P t; *- ma l 

value of S. is S v 
2 u2 

Inspection of (46) shows that it is 
independent of B provided that B^ < B < B 2 . When the value 
of B is reduced to B^, (46) reduces to (44); that is 



B 



b i r i 



< X 2 < T 



B 



b 2 r 2 



/ 
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* 

and, since S 2 must be B/b 2 for all B < B^ , this range 
must be valid for all such B. 

The analysis in section b above showed that 
(44) gave feasible x 2 values as long as B _< B Q . The analysis 
in this section then suggests that B^ must be less than or 
equal to Bq if (46) and (44) are to be feasible and identical 
at B = B^. Now the range of (46) is feasible if 



b 2 S u2 . S u2 

b l r l - r 2 



Collecting like terms results in 



b l r l 



b,S o < [ r =rx: ] T 

2 u2 _ b^r^ + ^2 r 2 



which reduces to B^ < Bq 



In summary, if B 2 = ^i s u i > 



B 



and 



B^ = t> 2 S u 2 < B < B 2 then x 2 is any value in the range given 
* * 

by (46) and = B/b-^ , S 2 = S u2 



When B £ B^ 



we again 



A A 

have x 2 given by (44) and = B/b^ , S 2 = B/b 2 . 

Figure 7 shows the range of optimal x 2 as 
a function of B for this case. It is immediately evident 
that x 2 = B Q/b^r^ is optimal regardless of the value of B 



(2) B 2 - b 2 S u2 , 



B 1 = b l S u l anc ^ B^ < B < B 2 



By similar arguments to those given above in 

* * 

S 2 = B/b 2 . The x 2 range 

from (33) reduces to 



section (1) we get = S ^ and 
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FIGURE 7: The range of optimal as a function of B 

for the case where B x < B^ and B-j < B£ • 
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; , b, s , 

-Hi < X, < T - J-2± 
r l - 2 - b 2 r 2 



( 47 ) 



when < B < B2 . When B <_ B^ , S^ = B/b^ , 

k * 

S2 = b A> 2 and x 2 can be obtained from ( 44 ) . Figure 7 also 
illustrates the range of optimal x 2 values for this case. 

( 3 ) B ± = b 1 S ul = b 2 S u 2 = B 2 

In this case, the interval B^ < B < B2 

* k 

does not exist and therefore we have S^ = S u ^ , S2 = S^ 
and the X2 interval is given by ( 37 ) when B^ = B2 < B < B^ . 

k 

When B _< B^ then = B/b^ , S2 = B/T^ and the X2 

interval is given by ( 44 ) . Figure 8 shows the optimal range 
of X2 as a function of B. 

Having now examined all possible cases it can be 
seen that there are three major possibilities: 



First: the problem is 'unconstrained if B B^ and 

partially unconstrained (the range of X2 is limited)' 

"k k 

if max [B2/B x ] < B < B^ . Solutions for and S 2 

are the economic order levels S . and S 

ul u 2 

Second: If the storage constraint is between min[B^,BQ] and 
max[B 9 ,B ] then the solutions depend on the relation- 
ship of the parameters to each other and to B. They 
can be classified according to two ranges for B: 
namely, Bq < B < B^ and B^ < B < B2 • 
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FIGURE 8: The range of optimal x 2 as a function of B 

for the case where B x < B 2 and B-j = B 2 . 



46 



* * 

Solutions for and S 2 are: 







1 


2 






* 


A 


* 


A 


1. 


S 1 


= S 1 ; 


S 2 = 


S 2 




* 




* 




2. 


S 1 


= B/b 1 ; 


S 2 = 


B 0 /b 




* 




* 




3. 


S 1 


= V b l ; 


S 2 


= B/b 




* 




* 




4. 


S 1 


= S ul '• 


S 2 = 


B/b 2 




* 




* 




5. 


S 1 


= B/b 1 ; 


S 2 = 


S u2 



The first three solutions correspond to the B n < B < B 

U X 

* 

case and x 2 can only have one value. The last two solutions 

* 

correspond to the B-^ < B < B 2 case and x 2 can take on any 
values within a computed interval. 

Third: In all cases where B £ min [B^,Bq ] the space 

constraint is so restrictive that only one item can 

* * 

be stored at a time. Hence = B/b^ and S 2 = B/b 2 , 

and the interval for x 2 expands as B decreases. 

The details of the optimal solutions of the two-item case are 
summarized by the flowchart presented in Figure 9. 

C. THE n-ITEM CASE 

The n-item case can be easily formulated using the same 
approach as for the two-item case. The general n constraint 
equations can be written for the ' times at which the items 
arrive using the max functions. 
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★ ★ ★ 

FIGURE 9: Flowchart for determining S-j , and for the 

two-item deterministic order level model. 
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1 




FIGURE 9: (continuation) 
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2 




FIGURE 9: (continuation) 
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